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The coupled cluster calculations with accounting for relativistic eflects to study spectroscopic and 
chemical properties of element 120 (E120) are performed. Similar calculations for Ba are also done 
and they are in a good agreement with the experimental data. Dissociation energies of diatomic 
X-H and X-Au molecules, where X=E120, Ba, are calculated; for E120 they are found to be 1.5 ~ 2 
times smaller than those for Ba. 



INTRODUCTION 

At present all the relatively long-lived superheavy el- 
ements (SHEs) up to the seventh period of the Periodic 
table were synthesized in FLNR JINR (see [H, Q and 
references) and their synthesis was mainly confirmed in 
LBNL ^] and GSI 4]. The new challenge for nuclear 
physics is the synthesis of the elements of the eight pe- 
riod. 

According to its position in the Periodic table, element 
120 (E120) is supposed to be the s element and homo- 
logue of Ba and Ra. The earlier obtained ground state 
configuration of E120 at the Dirac-Fock level, 7s^7p^8s^, 
is confirmed in recent calculation Q using the correla- 
tion potential method. Relativistic stabilization of the 
valence s orbital 6, 7] leads in general to higher exci- 
tation energies from this shell and, as a consequence, to 
weakening the bond strengths in its chemical compounds. 
Therefore, E120 can resemble the noble gas rather than 
to the alkaline earth element. Provided that E120 is syn- 
thesized, its properties can be investigated experimen- 
tally by the gas-phase chromatography method. Such 
scheme was sucessfully employed in FLNR JINR 8] and 
GSI 3 for E112 and El 14 where gold served as adsor- 
bent. Theoretical calculations are required for planning 
similar experiments with E120 on the gold surface. 

We have performed calculations of excitation energies 
of the E120 atom and its cation, E120'^, as well as of spec- 
troscopic properties of the E120H and E120Au molecules 
compared to the analogous systems with Ba, for which 
some experimental data are available. 



CALCULATION METHOD 

To calculate properties of E120, the generalized rel- 
ativistic effective core potential (GRECP) method ^ 
was employed. For relativistic correlation calculations, 
the fuUy-relativistic Fock-space coupled cluster code with 
single and double cluster amplitudes (FS-RCCSD) 0,[l3| 
was applied. To calculate corrections on enlargement of 
the basis set and higher cluster amplitudes, the CFOUR 
11 1 and MRCC H, 13] codes were used. 



GRECP generation and atomic calculations 

The GRECP for El 20 was generated in the frame- 
work of the present work. To check performance of 
the GRECP method for such a heavy element, relativis- 
tic correlation calculations on the E120 atom and its 
cations were carried out. We used four-component cal- 
culations with Dirac-Coulomb-Breit (DCB) Hamiltonian 
and Fermi nuclear charge distribution {A = 304) as the 
reference calculations. This Hamiltonian takes into ac- 
count the great bulk of the relativistic effects including 
the relativistic corrections to the Coulomb interaction 
between electrons. The 1— 5s, 2— 5p, 3—5d, and 4—5/ 
shells were frozen from the ground state E120^"'" calcu- 
lation. The 28 electrons (occupying the 6spd, 7spd, 8sp 
shells) were correlated in the [12,14,14,13,12,11] basis set 
of 6-17si/2, 6-19pi/2, 6-19p3/2, ... 6-I6/111/2 numer- 
ical spinors localized in the same radial space region as 
the 6spd, 7spd, 8sp spinors. The correlations were taken 
into account with the help of the FS-RCCSD method. 
The closed-shell ground state of E120^"'" was the refer- 
ence state and the Fock-space scheme was 
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with electrons added to the 8s 1/2, 8pi/2, 8753/2, 7^3/2, 
and 7dr,/2 spinors (the relativistic configurations corre- 
sponding to the nonrelativitic 8p^, 7d'^, and 8p^7o?^ ones 
were excluded from the model space). 

The leading configurations and terms for the lowest- 
lying states of the E120 atom and its cations are pre- 
sented in the first and second columns of table [H The 
results of the semiempirically-fitted CI/MBPT calcula- 
tions with Dirac-Coulomb Hamiltonian 0, |3] are com- 
pared with our FS-RCCSD results in the third and fourth 
columns. They are in a reasonable agreement. The abso- 
lute errors due to decrease in the number of the correlated 
electrons, the neglect of a finite nuclear size or Breit in- 
teractions are listed in the last four columns. One can 
see that a finite nuclear size should be taken into ac- 
count, whereas the correlations with the 6spd electrons 
and Breit interactions can be neglected for the accuracy 
within 2kcal/mol« 700cm^^ for one-electron excitations. 

We have generated the GRECP for only 10 explicitly 
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TABLE L Transition energies (TE) from DCB/FS-RCCSD 
calculations of the lowest-lying states of the E120 atom and its 
cations for 28 correlated electrons with Fermi nuclear model 
and Breit interactions taken into account. Absolute errors 
of reproducing the TE in the different versions of the four- 
component calculations. All values are in cm^^. 



Number of corr.el. 



28 20 10 



28 



28 



TE TE 



Abs. errors 



Leading 




Ref. 








Point without 


conf. 


Term 










nucl. 


Breit 


8s;/2 


(J=l/2) 














8Pl/2 


(J=l/2) 


24851 


24551 


55 


-372 


1299 


-51 


74/2 


(J=3/2) 




25355 


43 


-151 


1714 


135 


74/2 


(J=5/2) 




27529 


20 


-481 


1661 


152 


8P3/2 


(J=3/2) 


38057 37643 


23 


-401 


1632 


62 


£120^+ 


(J=0) 


89931 


89601 


36 


-620 


1592 


77 



8sv2 


(J- 


0) 














8si/28pi/2 


(J- 


0) 


16061 


15328 


58 


-290 


975 


-45 


8S1/28P1/2 


(J- 


1) 


17968 


17382 


50 


-290 


1015 


-30 


851/27^3/2 


(J- 


1) 


23066 


22337 


30 


-145 


1242 


101 


8Si/27d3/2 


(J- 


2) 


23231 


22494 


24 


-249 


1212 


101 


851/27^5/2 


(J- 


3) 


23827 


23377 


12 


-418 


1186 


108 


8S1/28P3/2 


(J- 


2) 


25457 


25308 


26 


-264 


1186 


42 


8^1/27^5/2 


(J- 


2) 


27477 


27652 


15 


-401 


1347 


100 


8S1/28P3/2 


(J- 


1) 


27685 


28304 


3 


-422 


1024 


39 


8sl/2 


(J- 


1/2) 


47296 


47633 


11 


-458 


906 


40 



treated electrons of E120 (i.e. with the 7spd, 8sp,6f, 5g 



GRECP components) following the scheme [15[ just 
to reduce unnecessary computational efforts at the 
stage of correlation molecular calculations (see below); 
similar GRECP version was generated earlier for Ba 
[l6j . The results for 10 correlated electrons (occupy- 
ing the 7spd,8sp shells) in the [9,11,11,9,8,7] basis set of 
7-15si/2, 7-17pi/2, 7-17p3/2, . . . 6-I2/111/2 numerical 
spinors localized in the same radial space region as the 
7spd, 8sp spinors are presented in table |lll Transition en- 
ergies from the DCB /FS-RCCSD calculations with Fermi 
nuclear model and the absolute errors of their reproduc- 
ing in the GRECP calculations are tabulated in the third 
and fourth columns. 

Different approximations to the "full" GRECP calcu- 
lations are considered in the last four columns. If one 
neglects the difference between the outercore (Jsp) and 
valence (8sp) GRECP components, two extreme GRECP 
versions can be derived with the conventional semi-local 
RECP operator: only valence or only outercore GRECP 
components acting on both the valence and outercore 
electrons. These cases are referred to as the valence or 
core GRECP versions, respectively One can see that 
the full and valence GRECP versions are suitable for 
the accuracy of 2 kcal/mol« 700 cm~^ whereas the core 
GRECP version is not. It should be noted that the er- 
rors of neglecting the innercore correlations with the 6spd 
shells and the errors of the GRECP approximation are 
partly compensating each other, thus, it additionally jus- 



tifies our choice of the 10-electron GRECP for the present 
molecular calculations. 

The scalar-relativistic (SR) calculations, i.e. without 
spin-orbit (SO) part of the valence GRECP operator, are 
presented in the seventh column. The SO contributions 
are large and should be taken into account. The scalar- 
relativistic SCF calculations followed by the FS-RCCSD 
calculations with the SO part of the valence GRECP op- 
erator are presented in the last column. These errors are 
comparable with the errors of the valence GRECP cal- 
culations in the fifth column. It should be emphasized 
that the GRECP calculations in the 4— 6-th columns were 
carried out with the SO part at both the SCF and FS- 
RCCSD stages. The computational versions used in the 
last two columns (the SR SCF calculation followed by 
the scalar-relativistic or fully-relativistic coupled cluster 
study) are also used in the molecular correlation calcula- 
tions discussed below. 



TABLE II: Transition energies (TE) from DCB/FS-RCCSD 
calculations of the lowest-lying states of the E120 atom and its 
cations for 10 correlated electrons with Fermi nuclear model 
and Breit interaction are taken into account. Absolute errors 
of reproducing the TE with different versions of the GRECP 
calculations. All values are in cm~^. 



TE 



GRECP abs.errors 















val. 


Leading 










val. 


SR-SCF 


conf. 


Term 


DCB full val. 


core 


SR SO-RCC 


8S1/2 


(J= 


.1/2) 










8PV2 


(J= 


=1/2) 24142 247 306 


-2761 


9255 


393 


74/2 


(J= 


=3/2) 25241 399 436 


-1766 


1316 


358 


74/2 


(J= 


=5/2) 27071 389 447 


-1765 


-514 


447 


8P3/2 


(J= 


=3/2) 37194 326 392 


-1937 


-3797 


257 




(J= 


=0) 


88907 446 554 


-1701 


-739 


417 


85^/2 


(J= 


=0) 










8sv28pi/2 


(J= 


=0) 


15012 266 310 


-2354 


6960 


-100 


8SV28P1/2 


(J= 




17064 251 295 


-2341 


4908 


-171 


851/27^3/2 


(J= 




22207 366 401 


-1632 


585 


316 


8Sj/27d3/2 


(J= 


--2) 


22259 387 428 


-1620 


533 


483 


8Si/27d5/2 


(J= 


=3) 


22968 386 436 


-1617 


-176 


484 


8SV28P3/2 


(J= 


=2) 


25009 253 297 


-1842 


-3037 


384 


8Sl/27rf5/2 


(J= 


--2) 


27271 344 394 


-1813 


-262 


291 


8S1/28P3/2 


(J= 


= 1) 


27834 307 351 


-1820 


-1171 


563 




(J= 


=1/2) 47120 333 385 


-1661 


-687 


381 



To compare with the E12G transition energies, the cor- 
responding experimental data for Ba from Ref. 17[ are 
listed in table IIIIl One can see that the barium excita- 
tion energies are smaller in general. It also indicates that 
E120 will possibly be more inert in general than Ba. 



Molecular calculations 

In two-component molecular relativistic calculations 
and high-level correlation treatment only relatively small 
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TABLE IIL The experimental transition energies (TE) from 
Ref. 17] for the lowest-lying states of the Ba atom and its 
cations. All values are in cm~^. 







Exper. 


Leading 




TE 


conf. 


Term 


Ref. [17] 


6si/2 


(J=l/2) 


— 5> 


1 

5^3/2 


(J=3/2) 


4874 


5'^5/2 


(J=5/2) 


5675 


6pi/2 


(J=l/2) 


20262 


6P3/2 


(J=3/2) 


21952 




(J=0) 
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6sv2 


(J=0) 


— >• 


6Sl/25'^3/2 

6Si/25d3/2 +6Si/25rf5/2 


(J=l) 
(J=2) 


9034 
9216 


6Sl/254/2 

6Sl/25'^5/2 +651/25^3/2 

6sv26p1/2 


(J=3) 
(J=2) 
(J=0) 


9597 
11395 
12266 


6SV26PV2 + 6si/26p3/2 

6SV26P3/2 


(J=l) 
(J=2) 


12637 
13515 


6SV26pJ/2 + 6S1/26P1/2 


(J=l) 


18060 




(J=l/2) 


42035 



basis sets can be used for diatomics like E120Au. At 
the same time, rather large basis sets can be employed 
in scalar-relativistic calculations of diatomic molecules. 
Therefore, the following scheme for the basis set gen- 
eration was used in this work: (i) For each atom 
(E120, Ba and Au), a large set of primitive Gaus- 
sian functions capable of describing wave-functions of 
the ground and excited states of the corresponding 
atoms was generated. These basis sets will be re- 
ferred as LBas below. LBas(E120) and LBas(Ba) consist 
of 15s—, 15p— , 8d—, 8/—, 6g—, 6/1— type functions, which 
shortly can be written as [15,15,8,8,6,6]. (ii) Then scalar 
relativistic CCSD calculation is performed with the large 
basis set for an atom and its compound (E120, E120- 
H, E120-Au, etc.). (iii) Generation of a compact basis 
set of contracted Gaussian functions was performed in a 
manner similar to that employed for generating atomic 
natural basis sets 
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the atomic blocks from the den- 
sity matrix calculated at stage (ii) were diagonalized to 
yield atomic natural-like basis set. The functions with 
the largest occupation numbers were selected from these 
natural basis functions. The results obtained with the 
given basis set approximately reproduce those with the 
large basis set. Besides, the functions required for accu- 
rate reproducing the essentially different radial parts of 
the 7pi/2 and 7p3/2 spinors, as well as the 8pi/2, 8^3/2, 
and 7c?5/2 spinors, have also been included to the 
new bases. These compact basis sets will be referred as 
CBas. 

Finally, the following scheme to evaluate ionization po- 
tentials and dissociation energies of molecules was em- 
ployed: (i) Calculation using two-component Fock-Space 
coupled cluster method with single and double ampli- 



tudes in the CBas basis set. (ii) Calculation of corrections 
on enlargement of the basis set and contribution of triple 
cluster amplitudes by the scalar-relativistic coupled clus- 
ter method with single, double and non-iterative triple 
cluster amplitudes, CCSD(T), using LBas. (iii) Calcula- 
tion of corrections on higher (iterative triple and non- 
iterative quadruple) cluster amplitudes by the scalar- 
relativistic coupled cluster method with single, double, 
triple and non-iterative quadruple cluster amplitudes, 
CCSDT(Q), using CBas. 

Equilibrium internuclear distances in the considered 
diatomic molecules were calculated using the scalar- 
relativistic CCSD(T) method and LBas. 



RESULTS AND DISCUSSIONS 
Ionization potentials 

Some properties of E120, of which the first (IPl) and 
second (IP 2) ionization potentials are examples, are con- 
sidered here in comparison with the corresponding prop- 
erties of Ba. As was described above, the two-component 
FS-RCCSD method was used to calculate the main con- 
tributions to IPl and IP 2. These calculations were 
performed in compact basis sets, CBas, consisting of 
5s, 6p, Ad, and 3/ functions for E120 and of 5s, 5p, Bd, and 
3/ for Ba. The computed values for E120 are IP1= 47236 
cm~i and IP2= 89061 cm-\ while for Ba IP1= 42340 
cm~^ and IP2= 80326 cm^^. Contributions from en- 
largement of the basis set up to 15s, 15p, 8d, 8/ and &h 
([15,15,8,8,6,6]) functions and non-iterative triple cluster 
amplitudes for E120 are -462 cm^^ for IPl and 31 cm~^ 
for IP2 [2^. For Ba these values are -579 cm~^and -65 
cm~^. Contributions of higher cluster amplitudes esti- 
mated using CCSDT(Q) method and [15,15,8,8,6] basis 
set {h functions were excluded) are negligible (less than 
30 cm~^). The final values of the ionization potentials 
for E120 and Ba are given in table IVTl together with the 
corresponding experimental values for Ba. 

The theoretical uncertainty of the ionization potentials 
of E120 can be estimated from the corresponding atomic 
calculations (tablesHlandUB and is suggested to be within 
1 kcal/mol«350 cm^^. 



X-H dissociation energies, X=E120, Ba 

In order to estimate the stability of compounds of E120 
compared to those of Ba we have first considered the dis- 
sociation energies of the corresponding hydrides since the 
experimental data are available for BaH [3]). Another 
often considered characteristic of SHE is the dissociation 
energy of its fluorides, X-F. However, the X-F bonding 
is not so illustrative qualitatively because almost all the 
elements (except light noble gases) are known to react in 
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a fluorine atmosphere yielding rather stable fluorides. At 
the same time the dimer systems such as Ba2, Hg2, Xe2, 
El 122 are all the van der Waals systems with small dis- 
sociation energies. By contrast, the ground state of the 
XeH molecule is not observed in the gas phase, whereas 



BaH was obtained and characterized 

To calculate E120H and BaH, the scheme similar to 
that for the calculation of the ionization potentials was 
used. Compact basis sets for E120, Ba and H were 
[5,6,4,2], [5,5,3,2] and [4,3,1], respectively. Large basis 
sets for E120 and Ba were [15,15,8,8,6], i.e. without h- 
functions. The aug-cc-pvqz [23| basis set was used as 
the large basis set for H. To exclude the basis set super- 
position errors, the diatomic molecules and atoms were 
calculated in the same two-center basis, i.e. the counter- 
poise corrections Q were used. 

The dissociation energy of E120H calculated within 
FS-RCCSD using the CBas(E120) and CBas(H) basis 
sets is 8258 cm~^, while the correction on the large basis 
and triple cluster amplitudes is -261 cm~^. For BaH the 
former contribution is 17061 cm^^, while the correction 
gives -430 cm~^. 

The calculated spectroscopic properties of the BaH 
molecule are in a good agreement with the experimen- 
tal data (see table llV| ). 



TABLE IV: Calculated and experimental harmonic frequency 
(we), vibrational anharmonicity {weXe), dissociation energy 
(Do) and equilibrium distance (Re) of the ^^^BaH molecule. 



Calculation Experiment 


We, cm~^ 


1158 


1168 [20,1 


WeXe, cm~^ 


14.1 


14.5 [2C[| 


Do, cm~^ 


16053 


15728 [19J 


R„ A 


2.238 


2.232 [20] 



The final values for the dissociation energies, D^, of 
E120H and BaH are given in table IVII The equilibrium 
internuclear distance in BaH is on 0.2 A shorter than that 
in E120H. 

The theoretical uncertainty of the dissociation energy 
of E120H is estimated to be 500 cm~^. 

It follows from table Ell that the E120-H bond is sig- 
nificantly weaker than that in BaH. Partly it can be ex- 
plained by the presence of lower-lying excited states in 
the case of Ba and its cation (see table IIIll) . Contribution 
of these states to the BaH chemical bond leads to its ad- 
ditional stabilization. In turn, the excited states in the 
case of E120 are lying significantly higher (see table HI]); 

To check this viewpoint two series of scalar-relativistic 
calculations of the BaH and E120H dissociation energies 
were performed: (i) with d~ and /—type basis functions 
in basis set (CBas); (ii) without d and / functions. The 
latter calculation prevents participating the d orbitals of 
Ba and E120 in chemical bonding of their hydrides. Note, 
however, that it is just test calculation because exclu- 
sion of d basis functions also prevents, e.g., correlation of 



5p— electrons of Ba and 7p— electrons of E120 into these 
states. The results are given in table IVl 



TABLE V: Calculated dissociation energy of BaH and E120H 
using scalar-relativistic Hartree-Fock (HF) of CCSD(T) levels 
of theory in CBas basis set and CBas with excluded d— , /— 
functions (CBasNoD) 



Basis 


CBas 


CBasNoD 




Method 


HF 


CCSD(T) HF 


CCSD(T) 


De(BaH), cm"^ 


13605 


15957 7877 


10677 


De(E120H), cm" 


2864 


6865 -68 


5264 



It follows from table |V] that d functions (and higher 
harmonics) significantly contribute to bonding of the 
monohydrides under consideration. At the CCSD(T) 
level of theory, this contribution is 5280 cm~^ to the 
Ba-H bond, while is only 1601 cm^^ in E120H. These 
values confirm the qualitative discussion above based on 
the atomic transitions. However, even without d basis 
functions the BaH bond energy at the CCSD(T) level is 
twice stronger than that in E120H. This observation is 
in a qualitative agreement with the fact that the states 
with the valence sp (p) configurations in E120 (E120+) 
lies higher than the corresponding states in Ba [Ba^). It 
should be noted that the bonding in E120H is in essence 
due to the correlation effects when the CBasNoD basis 
set is used. 



X-Au dissociation energies, X=E120, Ba 

As the first stage of modelling interaction of E120 ele- 
ment with gold surface, the simplest comparative model, 
E120Au vs. BaAu, is considered here. 

The 19-electron GRECP was used for Au. Thus, 29 
electrons were treated in the correlation calculation. For 
E120Au, the relativistic two-component FS-RCCSD cal- 
culation in the Au[7,7,4,2] and E120[5,5,2,l] basis sets 
gives the dissociation energy of 11732 cm~^. Correc- 
tion on the larger basis set ([15,15,8,8,7] for Au and 
[15,15,8,4] for E120) is 510 cm~^, contribution of non- 
iterative triple cluster amplitudes [25] is 2264 cm"-'^ and 
correction on higher amplitudes (calculated as the differ- 
ence between CCSDT(Q) and CCSD(T) energies in the 
compact Au[6,6,4,2,l] and E120[6,5,2,l] basis sets) is less 
than 100 cm~^. Similar calculations were performed for 
BaAu. The final calculated values are given in table IVII 
The equilibrium distance for BaAu is on 0.1 A shorter 
than that in E120Au. 

The theoretical uncertainty of the dissociation energy 
of E120Au is estimated to be 1000 cm~^ 
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TABLE VI: Calculated properties of E120 in comparison with 
Ba: the first and second ionization potentials (IPl and IP2), 
dissociation energies (-De) of hydrides and aurides (X-Au) and 
the corresponding equilibrium distances (Re)- 



E120-calc Ba-calc Ba-exp 


IPl, cm"' 


47046 


41932 


42035 [24] 


IP2, cm"^ 


89286 


80442 


80686 [24] 


De(X-H), cm~" 


7997 


16632 


16308[19] 


i?e(X-H), A 


2.45 


2.24 


2.23 [20| 


De(X-Au), cm~^ 


14428 


22608 




Re{X-An), A 


3.1 


3.0 





Conclusion 

Properties of E120 and its compounds are considered 
in comparison with their Ba analogues. The monohy- 
dride and monoauride of E120 are found to be less stable 
than the corresponding analogues of the Ba compounds. 
Nevertheless, E120 can be rather considered as a "typi- 
cal" representative of the second group. 
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